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Abstract 

Pirst we consider a unidirectional flux u) of vehicles each of which 
is characterized by its 'natural' velocity v drawn from a distribution 
P{v). The traffic flow is modeled as a collection of straight 'world 
lines' in the time-space plane, with overtaking events represented by a 
fixed queuing time r imposed on the overtaking vehicle. This geomet- 
rical model exhibits platoon formation and allows, among many other 
things, for the calculation of the effective average velocity w = (j){v) 
of a vehicle of natural velocity v. Secondly, we extend the model to 
two opposite lanes, A and B. We argue that the queuing time r in 
one lane is determined by the traffic density in the opposite lane. On 
the basis of reasonable additional assumptions we establish a set of 
equations that couple the two lanes and can be solved numerically. It 
appears that above a critical value Wc of the control parameter oj the 
symmetry between the lanes is spontaneously broken: there is a slow 
lane where long platoons form behind the slowest vehicles, and a fast 
lane where overtaking is easy due to the wide spacing between the 
platoons in the opposite direction. A variant of the model is studied 
in which the spatial vehicle density p rather than the flux Co is the 
control parameter. Unequal fluxes u)a and in the two lanes are 
also considered. The symmetry breaking phenomenon exhibited by 
this model, even though no doubt hard to observe in pure form in 
real-life traffic, nevertheless indicates a tendency of such traffic. 

Keywords: bidirectional traffic, overtaking traffic, geometric model, 
spontaneous symmetry breaking 
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1 Introduction 



Traffic issues arise in many physical systems, ranging from intra- 
cellular traffic to road traffic pQ. The development of simple traffic 
models, such as exclusion processes, has led to a deeper understand- 
ing of how such out-of-equilibrium systems are driven (for example how the 
introduction of a reaction time can induce metastability [2]). Much effort has 
been devoted to unidirectional flows, sometimes on several parallel tracks (or 
'lanes' in the road traffic vocabulary) [21 IH El El El El IS]- However, many 
applications involve bidirectional flows. It is only recently that models for 
bidirectional traffic have been proposed. 

1.1 Bidirectional traffic and symmetry breaking 

In one set of models the two lanes are represented by a single chain of lat- 
tice sites. Oppositely moving entities (vehicles, particles, molecular motors, 
...) share the same lane and certain exchange rules have to be defined for 
when they meet. Recently such a model has been solved exactly for a whole 
distribution of hopping rates [10] . In a variant called the "bridge model" , op- 
positely moving particles have to slow down when they meet, as on a narrow 
bridge [m Ha [131 HI ttSl [16]. 

In another set of models the two lanes are represented explicitly. When 
both lanes are accessible to particles moving in either direction, jams may be 
formed [T71 [TH [191 [20]. By contrast, the traffic flow is quite efficient when 
oppositely moving particles are confined to different lanes [2T]. Even in that 
case, however, residual interactions between particles on different lanes can 
be taken into account (e.g. [22]). For example, in p3j a particle (or vehicle) 
has to slow down when it crosses another one on the other lane; however, 
the authors restrict themselves to the special case where one of the lanes 
contains only one vehicle. The influence on two-lane bidirectional traffic of 
quenched disorder (in the sense of statistical physics) [21] , or of the dynamics 
of the track itself (in the context of intracellular traffic) [25], has also been 
considered. All the models that we have mentioned so far are defined on a 
discrete space and the motion consists of hopping events. 

One phenomenon of interest in traffic models is the spontaneous breaking 
of symmetry, whether it be between two - in principle equivalent - lanes or 
between two flow directions. Symmetry breaking in discrete traffic systems 
has been reported for two-lane unidirectional traffic [261 123 [28] and for one- 
lane bidirectional traffic (e.g. [21], or in the bridge model [HI [121 [131 El 
[T5| [IS]). In [2^1 [3D1 IS] symmetry breaking occurs for two-lane bidirectional 
traffic, but the only interaction between the lanes is through narrow road 
entrances - which may be seen as localized bridges - at the two road ends. 
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Figure 1: A two-lane road has vehicles (black circles) on lane A moving to the 
right and vehicles (red circles) on lane B moving to the left. Arrows of variable 
length correspond to a distribution of velocities. Strings of circles joined together 
indicate platoons of faster vehicles queuing behind a slower one. In the center of 
the figure a black vehicle (filled black dot) is overtaking a platoon moving in the 
same direction. Platoon formation, overtaking events, and interaction between the 
two lanes are the phenomena that we seek to describe in this work. In the actual 
model the vehicles and the platoons are point-like. 

In the present work our motivation is to construct a description of the 
traffic situation depicted in figure [1] using a minimum number of mathemat- 
ical ingredients. This leads us to propose a new two-lane model which is 
continuous in space and time. It is simple enough to allow for theoretical 
analysis, yet retains the essentials of real-life traffic for low to moderate ve- 
hicle densities. A continuous distribution of velocities is considered. On a 
given lane vehicles with different velocities may overtake each other provided 
there is enough space on the other lane. 

The paper consists of two parts. In the first part, section [2], we consider 
only a single lane of the two-lane road and the infiuence of the other lane 
comes in only through a postulated input function. In the second part, 
section [3l we consider two opposite lanes coupled together and determine 
the input function. We report spontaneous symmetry breaking between the 
lanes and are able to demonstrate analytically how it results from the bulk 
evolution rules. In sectionlHwe show numerical results. In section[5]we collect 
a number of comments on the exact mathematical status of the present one- 
and two-lane models. Section [6] is our conclusion. 

1.2 The present model 

1.2.1 First part: single-lane model 

The single-lane model of section [2] is purely geometric. It is governed by 
the rule (see figure [2]) that a vehicle with velocity v can overtake a slower 
one with velocity v' < v only after it has stayed queuing behind that slower 
vehicle during a queuing time time t{v'). Here t{v') is the waiting time 
that elapses, on average, before the queuing vehicle finds a sufficiently long 
free time interval to execute the overtaking maneuver. In the terminology of 
physics it is an interaction constant between the vehicles that depends only 
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Figure 2: The basic overtaking event of a vehicle of velocity v' by a vehicle of 
higher velocity v. During the time t{v') necessary for this event the two vehicles 
form a 2-platoon. As a result the trajectory of the overtaking vehicle is displaced 
to the right by an amount tdeiiv,v') = (1 — v' /v)t{v'). 

on the velocity of the slower vehicle. A special case is obtained by taking 
the queuing time constant, t{v) = f. In the single-lane model t{v') is an 
arbitrarily postulated function. Its expression should be derived, in principle, 
from the traffic density in the opposite lane, which however at the level of 
the single-lane model is not described explicitly. 

During the time interval t{v') the vehicle of velocity v' is the leader and 
the vehicle of velocity v is the follower. After the time interval t{v') the faster 
vehicle resumes its original, or 'natural' velocity v. By a simple geometric 
argument in the time-space plane, the vehicle with the higher velocity incurs 
a time delay (a displacement to the right of its trajectory in the time-space 



At high traffic densities more than one vehicle may be queuing simultane- 
ously behind the same slower vehicle and hence overtaking events will occur 
that involve three or more vehicles. The trajectories for such many-vehicle 
events may be constructed with the aid of an appropriate interpretation of 
the same queuing time rule (see section [2]). The rules then allow for the 
construction of the full set of individual vehicle trajectories in the time-space 
{tx) plane. We call this model 'geometric' because the resulting vehicle tra- 
jectories form a pattern of straight lines (a sample pattern will be shown in 

^In functions with two velocity arguments, the first argument (velocity of the vehicle 
that overtakes) will always be larger than the second one (velocity of the vehicle being 
overtaken) . 



plane) equal tcEl 




(1.1) 



4 



section [21 figure 5). 

The construction shows that at any point of time each vehicle is either 
a follower (that is, it is being delayed behind another vehicle) or, if not, is 
a leader and advances at its natural velocity accompanied by any number 
(possibly zero) of followers. A leader accompanied by n — 1 followers will be 
called an n-platoon; a 1-platoon, therefore, is a single free vehicle. 

In space-time diagrams like those of figure |2] a moving platoon is repre- 
sented by coinciding straight-line segments. In our model both single vehicles 
and platoons of vehicles are point-like: they have zero length in space. As a 
consequence, the model can be expected to apply to low and moderate ve- 
hicle densities p{v), but not to very high ones. For this single-lane problem 
we will determine, in particular, the two quantities most characteristic of a 
stationary traffic flow. These are the spatial density pica.d{v) of leaders (which 
is also the platoon density) and the spatial density Pfoii(f ) of followers, for 
a given velocity v. The expressions for these quantities will depend on the 
queuing time function t{v) and on the boundary conditions. The latter may 
be imposed either at point of entry x = of the road section (for a road 
section with open-ended boundary conditions), or at an initial time t = 
(for a "ring": a road section with periodic boundary conditions); the two 
possibilities give rise to slightly different mathematics. 

A few words are in place here about other geometric single-lane models 
occurring in the literature. As an example we cite Ben-Naim, Krapivsky, 
and Redner [52] , in whose model a vehicle or a platoon that catches up with 
another one assumes that one's velocity without ever overtaking it; the model 
therefore describes clustering into platoons that can only grow. It has no 
stationary state but is appropriate for the analysis of platoon formation, for 
which the authors derive scaling laws. In an extension of the same model Ben- 
Naim and Krapivsky [321 133] consider overtaking based on vehicles having 
an "escape rate" from the platoons; there is then a stationary state which 
may be analyzed by a Boltzmann type approach. 

1.2.2 Second part: two-lane model 

In section [3] of this paper we consider two opposite lanes called A and B, 
as represented in figure [1], each one described individually by the single-lane 
model of section [2] with queuing times r^(f) and tb{v). The queuing time 
ta{v) is supposed to have its origin in the traffic density in lane B and vice 
versa. On the basis of additional hypotheses we express ta{v) explicitly 
in terms of the traffic density in lane B, and vice versa. This essentially 
amounts, in physical parlance, to establishing a mean-field type interaction 
between the traffic flows on the two lanes. It is again possible to study the 
resulting model analytically. The simplest case occurs when there is full sym- 
metry between the lanes A and B. Let a parameter oj represent the intensity 
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of the traffic flow. The most interesting result is that in that case, when u 
crosses a critical value Wc, the symmetry between the lanes is spontaneously 
broken: one lane, say A, will have a high density of short platoons with 
on average a higher velocity than the second lane, which will have a lower 
density of longer platoons. The platoon density difference pA.iead — Ps.iead 
plays the role of the order parameter, as will be seen in figure [HI We de- 
scribe this phase transition and extend our analysis and simulations to the 
non-symmetric case. 

Because of the nonlocal (mean-field) character of the lane-lane coupling, 
this two-lane model, contrary to the single-lane models on which it is based, 
does not allow any longer for a geometric representation. That is, there is 
no corresponding model of individually interacting vehicles. In a companion 
paper [35] we will present a truly microscopic geometric two-lane problem 
and compare its simulation results to the theory of this work. 

2 A single-lane model 

We will ffist complete the definition of the single-lane model described in 
section 11.21 This model has unidirectional traffic and a basic overtaking 
event as pictured in figure |21 The delay time r(f ) is a fixed given function. 
We will associate with each natural velocity v an effective or average velocity 
w = 0(f), then find an equation for the function (p, and show how the platoon 
statistics may be derived from its solution. 

2.1 Preliminaries 
2.1.1 Many- vehicle events 

We must begin by briefly dwelling on overtaking events involving three or 
more vehicles. For these we will adopt the rule that queuing times, and 
hence incurred delays, are additive. Figure |3] shows an example of a three- 
vehicle event. The vehicle of velocity v overtakes successively two slower 
vehicles of velocities Vi and V2 and as a result incurs a time delay = 
tdeiiv, Vi) + tdei{v, V2) with tjei givcu by (11. ip . One may note that the vehicle 
of velocity f 2 also overtakes the one of velocity vi and that there is a time 
interval during which the latter has two followers (there is a 3-platoon). 

In figure H] a slightly more complicated three- vehicle overtaking event is 
depicted. The vehicle of velocity v spends two disjoint time intervals and 
Tfe queuing behind the one of velocity V2, but such that Ta + = ^(^2), thus 
satisfying the rule imposing the total queuing time. 

One may construct arbitrary many-vehicle events respecting the same 
rule by drawing the trajectories of vehicles of successively higher velocities. 
As a consequence, in an event where a vehicle of velocity v overtakes k other 
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Figure 3: Example of a three- vehicle overtaking event with the time and space 
axes as in figure [2j First, the vehicle of velocity vi is overtaken by those of velocities 
V2 and V during partially overlapping time intervals both of length t{vi). During 
the overlapping part of those intervals there is a 3-platoon. Next, in a separate 
event of duration t{v2), the vehicle of velocity V2 is overtaken by the one of velocity 
i;. As a result the latter has incurred a time delay tdei(f^)^^i) + idei(f^; ^^2)- 

vehicles, its outgoing trajectory in the tx plane will again be a straight line 
displaced by an amount t^"* = J2i=i^dei{v,Vi) to the right with respect to its 
incoming trajectory. This is exhibited in figure [5l 

In figure |6l finally, we show a sample set of trajectories on a given lane 
section during a given time interval. 

2.1.2 Boundary conditions 

In order for the single-lane problem to be well defined it has to be supple- 
mented with boundary conditions. In the tx plane these may be applied 
either along the space axis or along the time axis. Both possibilities have 
merits, as we will briefly discuss now. 
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Figure 4: A different example of a three- veiiicle overtaking event; the time and 
space axes are as in figure [2j When the vehicle of velocity v has queued for a time 
Ta behind the one of velocity V2, the two catch up with the still slower vehicle of 
velocity vi. During a time interval t{vi) the two simultaneously spend the required 
queuing time t{vi) behind that vehicle, which they then simultaneously overtake. 
Following that, the vehicle of velocity v spends the remaining time = t[v2) — Ta 
queuing behind the one of velocity V2, after which it overtakes it. The net result is 
that the vehicle of velocity v has incurred the same time delay tdei(^^) vi)+t(i^\{v^ V2) 
as in figure [3l 
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Figure 5: A vehicle of velocity v overtakes k slower vehicles. Its final trajectory 
suffers a time delay which is the sum of the time delays incurred due to k 
individual overtaking events. The full sequence of partially overlapping ovcrtakings 
is a complicated event that occurs inside the dashed disk. It has not been presented 
explicitly: only the final outcome is relevant in our description. 
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Figure 6: Sample pattern of vehicle trajectories in the tx plane. The heavy line is 

the trajectory of a vehicle of natural velocity v. When averaged over a long time 
interval, it appears to have an effective velocity equal to the slope of the dotted 
line and denoted hy w = 4>{v). 
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Boundary condition in x = for all times (open-ended boundary condi- 
tions). - One may consider an open-ended lane section of coordinates < 
X < L and impose, for each f , the rate u{v)dv at which vehicles with natural 
velocity between v and v -\- dv enter the lane section at the point x = 0. We 
will set 



where now P{v)dv is the probability that an entering vehicle picked at ran- 
dom has a velocity in [v, v + dv]. In x = L the vehicles leave the lane section 
freely. We will refer to this situation as having open-ended boundary con- 
ditions (OBC). Because of the continuity equation, in a stationary state all 
vehicles that enter the lane at x = at some rate must also cross any other 
observation point a; > at the same rate. Hence the control function u{v), 
once given for x = 0, will be the same at any other point in space. 

Boundary conditions at t = in all of space (periodic boundary condi- 
tions). - One may consider, especially in the context of simulation [35] or 
of theoretical work, a circular lane of circumference L. It is then natural to 
impose at an initial time t = the spatial density p{v)dv of vehicles with 
velocity between v and v -\- dv. We will set 



where now R{v)dv is the probability that a vehicle picked at random at the 
initial time has a velocity in [v,v + df]. We will refer to this situation as 
having periodic boundary conditions (PEG). Since vehicles cannot leave or 
enter the ring, the control function p{v), when given at some initial time 
t = 0, will be the same at all later times t > 0. 

We will refer to the probability densities P{v) and R{v) as the 'control 
functions', and to the parameters u and p that multiply them as the 'con- 
trol parameters'. As shown above, both control functions are actually bulk 
quantities. 

The boundary conditions as we have formulated them, whether at x = 
or at t = 0, do not specify any correlations, or the absence thereof, that 
may exist between the positions and/or the velocities of the vehicles on these 
boundaries. We will briefly return to this point in section [51 In any case, we 
will make the 'physically reasonable' assumption that there is a penetration 
length Xpen in the case of OBC (and a relaxation time t^ei in the case of 
PEG) such that for x ^ Xpen (for t ^ trei) a state is reached that is both 
homogeneous in space and stationary in time. We will not attempt here to 
express Xpen in terms of u{v) and r(f), or trei in terms of p{v) and t{v). 





(2.1) 





(2.2) 
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2.1.3 Natural velocity v and effective velocity w 

The first question that comes to mind is suggested by figure [61 Suppose the 
traffic is in a stationary state. What will then be the effective velocity that 
is, the velocity averaged over a sufficiently long time, of a vehicle that has v 
as its natural velocity? Let us assume that there exists a relation yielding w 
clS cl function of v and write it as 

w = 0(f), V = ip{w), (2.3) 

where i/j is the inverse of 0. We would like to determine (j){v) for all v > vo- 
Let us denote the minimum velocity (which may be zero) by vq , that is, 
P{v) = (for OBC) or R{v) = (for PBC) when v < vq. The vehicles with 
V = Vq follow unperturbed trajectories: since they cannot overtake any other 
vehicle, they have straight world lines in the tx plane. We therefore have 
Wq = 0(fo) = Vq. Furthermore 0(f) should increase with v. 

Below we will find the expression for 0(f) and see that it is boundary 
condition independent. In subsection 12.21 we will discuss the case of OBC; in 
section [23] we will then describe succinctly what changes in the case of PBC. 



2.2 Open-ended traffic lane 
2.2.1 General 

Without the delays due to overtakings (that is: for r(f) = 0, or equiva- 
lently: in the absence of 'interaction' between the vehicles) the spatial density 
p^'^''(f)df of vehicles of natural velocity in [f,f + df] is given by 

p(0)(^)d^;= ^df. (2.4) 
V 

The factor here indicates that a fast vehicle spends less time in the 
lane section and therefore contributes less to the density. The total vehicle 
density, again for r(f ) = 0, is therefore 

p(o)= rdfp(°)(f)= Tdf^. (2.5) 

J vo J vo ^ 

This work is concerned with the modifications that occur due to the 
interaction between the vehicles. 

In the stationary state, at any fixed observation point vehicles of effective 
velocity w will pass with a probability of a)(w) per unit of time. We may find 
Cj{w) by first noting that 

Ci;(f )df = Cj{w)dw, (2.6) 
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which is the continuity equation. In fl2.6p and henceforth, tildes will indicate 
quantities expressed as functions of the effective velocity w. It follows from 
(12:61) and Q that 

u{w) = uj{i){w))ip'{w). (2.7) 
The spatial density p{w) of vehicles with effective velocity in [w, w + dw] is 

piw)dw = ^dw = ^dv. (2.8) 

W (j)[v) 

The spatial density p{y)dv of vehicles of natural velocity in [t>,f + dv] is 
defined by 

p{v)dv = p{w)dw. (2.9) 

We recall that each vehicle has a natural velocity that is fixed once and for 
all. The density p{v) takes into account all vehicles of natural velocity v 
regardless of their actual velocity at the moment. It thus appears from (12.81) 
and O that 

= (2.10) 

(j){v) 

This is a key equation of the present model, valid in its stationary state. 

By integrating (I2.10p over v we see that the total vehicle density in the 
interacting system is equal to 



Vo 



At this point 0(f) is still unknown. 



2.2.2 Deriving an equation for (j){v) 

We will now derive an equation from which 0(f) can be solved. The basic idea 
is provided by the fact that the trajectory of a vehicle with a given velocity is 
affected ('renormalized') only by those of vehicles with lower velocities. This 
idea can be exploited as follows. 

We continue to suppose the system is in a stationary state and consider it 
during a time interval [0, T]. In the limit of large T all vehicles will have well- 
defined effective velocities. Consider a marked vehicle of natural velocity v, 
and hence of effective velocity w = 0(f), starting in the origin x = at time 
t = 0. During the interval [0,T] any vehicle with effective velocity w' < w 
initially present in the spatial domain [0, {w — w')T] will be overtaken by the 
marked vehicle. Let N{w,w')dw' be the number of vehicles with velocity in 
[w',w' + dw'] overtaken by the marked vehicle (see footnote ^). In the limit 
of large T, again, this number will be given by 

N{w, w')dw' = {w- w')T X p{w')dw'. (2. 12) 
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We will write f{w) = t{iIj{w)) for the queuing time necessary for overtaking 
a vehicle of velocity w. The total duration of the time intervals occupied by 
the overtaking events (12.121) amounts to a time Tfo\i{w,w')dw' given by 

fMiiw,w')dw' = f{w')N{w,w')dw'. (2.13) 

During this time the marked vehicle is a follower traveling only at velocity 
ipiw') and hence covers a distance dfon{w,w')dw' given by 

dioiiiw, w')dw' = f{w')N{w, w')ij{w')dw'. (2.14) 

By integrating this last expression over w' and using (12.121) we obtain the sum 
dfoiiiw) of the distances traveled by the marked vehicle during its overtaking 
events, 

di^ii{w) =T dw'{w- w') p{w')'ip{w')f{w'). (2.15) 

JWQ 

The remaining time T\ea.d, not used for overtaking events, is 

fyeUw) = T-T dw'iw-w')piw')fiw'). (2.16) 

J wo 

During the time Tiead(w) the marked vehicle drove at its natural velocity v 
and covered a distance c?iead(w^) = vTiei^^^w) , where our notation is hybrid. 
We now obtain the effective velocity w of the marked vehicle as 

W = (dleadiw) + Jfoll(w))/T 

= v — v dw' {w — w') p{w') f{w') + / dw' {w — w') p{w')^jj{w') f{w'). 

J wo J wo 

(2.17) 

Employing (12. 3p . (12.101) . and (12. 7p in (12.171) and eliminating the variables w 
and w' , we reexpress (I2.17P as 



0(f) = V — dv' uj{v') t{v'){v — v') 

J ' 

which may be rewritten as 



Do 



0(f') 



(2.18) 



f + V dv' ijj{v')t(v'\(v — v'^ 
^ 1 + £dv' u{v')t{v'){v - v')[(f){v')]-^ ' ^' ' 

This is the equation for 0(t> ) that we sought. It depends on the shapes of 
the velocity distribution u){v) and the queuing time distribution t{v). 
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2.2.3 Solving the equation for 0(f) 



The equation for 0(f) may be readily solved. The easiest way to do this is 
to observe that it may be written as a linear equation in terms of l/0(f ). In 
appendix El it is shown that 



One easily checks that 0(fo) = fo as had to be expected. One also verifies 
that r(f") = leads to 0(f) = f , as had to be the case, and that r(f") = oo 
gives 0(f) = fo, which is intuitively obvious. 

There are a few cases in which the integrals in fl2.20p can be carried out 
exactly. One of them is the case where r(f) = f is a constant and w(f) is a 
block distribution, constant for f o < f < f m and zero elsewhere. We present 
it as an example in appendix [HI However, we will not pursue such special 
cases here. In the numerical section H] we present a figure showing the typical 
behavior of 0(f) for two different traffic ffuxes. 

In cases where the velocity distribution uj{v) is a sum of delta peaks, the 
reasoning performed above leads to discretized versions of equations fl2.19p . 
fl2.20p . and other results below. We brieffy present these in appendix O 

2.2.4 Platoons, leaders, and followers 

An n-platoon consists of n — 1 vehicles all blocked behind a single 'leader' 
advancing at its natural velocity, say v'. The n—1 'followers', although having 
natural velocities larger than v', are also advancing at the velocity v' of the 
leader. Each follower is in the queuing interval r after which it will overtake 
the leader. It will be convenient to refer to n as the 'platoon length', even if 
in the present model platoons are point-like. A freely advancing single vehicle 
without followers will be called a '1-platoon'. In this subsection we study the 
spatial density of platoons characterized by a given velocity and/or platoon 
length. Our approach will be to first study the statistics of the followers and 
then derive from it the statistics of the leaders. 

Followers. - Let /foii(w) be the fraction of its time that a vehicle of 
natural velocity v = tp{w) spends following slower vehicles. From equation 
f l2.16p we see that 




n -2 



(2.20) 





(2.21) 
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Similar reasoning, but without integrating on f ', shows that (see footnote 



fiv,v')dv' 



- 1 



uj{v')r{v')dv' , V > v', 



(2.22) 



is the fraction of its time that a vehicle of natural velocity v spends following 
vehicles with velocities in [v',v' + df']. We will now pass from these time 
fractions to spatial densities. 

Equation f l2.10p gives the density p{v) of vehicles of natural velocity v. It 
follows that in a stationary state the density of following vehicles of natural 
velocity v, to be called Pioii{v), is given by 



Pfoii(^^) = p{v)ffon{v). 
Hence, using f l230|) and (jMI]) in fl2:23|l we obtain 



PM\[V) 







/ dv' 









uj{v')t{v'). 



(2.23) 



(2.24) 



Let pfou denote the total density of following vehicles. This quantity is ob- 
tained as the integral on v of equation (]2.24p . Hence 



PfoU 



df 



V<P{v') 



uj{v')t{v'). 



(2.25) 



Let us define, for all v > v' 



p{v, v') dv dv' = spatial density of vehicles of natural velocity in [v, v + dv] 
following a leader of velocity in [v', v' + dv']. (2.26) 

By interpreting (12.251) we see that 

p{v,v') = ^,^,r{v') [<p{v) - <P{v')], V > v'. (2.27) 

0(f) (j){v') 

The two-velocity density function p{v, v') together with the expression for 
(j){v) are two fundamental quantities in the present model. We have the 
obvious relation 

Pfoii(tO = / dv'piv^v'). (2.28) 



In addition to f l2.28p there is a second integral that we can construct by 
projecting p(t>, v') onto the v' axis. 



dv p{v, v'). 



(2.29) 
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It has the interpretation of the total density of vehicles that follow a leader 
of natural velocity v' . We observe that 



PfoU 



VO J VO 



(2.30) 



which we record for later use. 



Leaders. - The quantity fxeadiy) = 1 — fiowiy) is the fraction of its 
time that a vehicle of natural velocity v spends not queuing: it may either 
constitute a 1-platoon all by itself or be at the head of an n-platoon. Upon 
setting correspondingly 



Plead (^^) = piv 

we find with the aid of (12.241) that 

Plead (^^) = P(f)/lead(^^) 



(2.31) 



1 



0(^7 ^ Jvo 

By integrating fl2.32p on v one obtains 



ct>{v'] 



- 1 



uj{v')t{v') 



Plead 



I VO 



P{v')r{v') 



We obviously have the relation 

Plead = P — Pfoll ■ 



(2.32) 



(2.33) 



(2.34) 



The quantities piead('y) and piead are also equal to the density of platoons of 
velocity v and the total platoon density, respectively. 

We may eliminate 0(f) from (I2.32p by substituting for it the solution 
f l2.20p . The resulting expression can be simplified importantly. In appendix 
Id] we show that after substantial rewriting it takes the form 



Plead l^^j 



u:{v) 



(2.35) 



which relates piead('^^) directly to u^v), without the intervention of 0(f). We 
retain fl2.35p for later use. It is employed, moreover, in appendix |B] to ob- 
tain the analytic expression for piead('y) in the case where u}{v) is a block 
distribution. 

Platoon length distribution. - We consider a leading vehicle having a 
velocity that we now call v (it was called v' above) and ask about the distri- 
bution of the number n of vehicles in its platoon. 
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At any instant of time, in a spatial interval of length L the total average 
number of leading vehicles of natural velocity in [v,v + dv] is Lpiead('y)dw = 
L[p{v) — pfoii(f )]dt'. The average total number of vehicles in this interval 
that are queuing behind leaders of natural velocity in [v, v + dv] is Lafon{v)dv. 
Hence a leading vehicle of natural velocity v is followed by an average number 
of vehicles i>{v) given by the ratio of these two numbers, 

p{v) - pioii{v) Plcad{v) 

The number 

n{v) = 1 + iy{v) (2.37) 
is the average length of a platoon of velocity v. We have the sum rules 

dv pic^d{v) i^{v) = pfoii, 

dv pic^div) n{v) = p. (2.38) 



Do 
OO 



Vo 



The platoon length n averaged on platoons of arbitrary velocities present in 
a given interval is equal to 

n = 1 + u 

Plei,d{v)u{v) 



OO 



1 + / dv' 

' VQ Plead 



Plead 



(2.39) 



where we used ([23S1), (I23DD, and IKM^ . 

A more detailed analysis of the composition of platoons is possible along 
the same lines. The total number of vehicles of natural velocity in [v,v + 
dv] queuing behind leaders of natural velocity in [v', v' + dv'] is equal to 
Lp{v, v')dvdv'. Hence a leading vehicle of natural velocity v' is followed by 
an average number i/(f , v')dv of vehicles of natural velocity in [v, v + dv] given 
by 

iy{v, v')dv = ,^ \ . ■ (2-40) 

p{v') - ploll{v') 

We have 

/■OO 

iy{v') = / dvu{v,v'). (2.41) 



Certain questions that one may ask require knowledge not only of the 
the average platoon length (12.411) . but of their probability distribution. Let 
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p„(f ) be the probability that a randomly picked leading vehicle of velocity v 
is heading an n-platoon (n = 1, 2, . . .). The work done so far does not allow 
us to write down an expression for Pn{v). If however we make the additional 
assumptioiJl that the followers are distributed randomly over the leaders, 
then Pn{v) is given by the Poisson distribution 

p„(^) = e--('')M^jll^, n = l,2,... (2.42) 

This distribution is required, for example, when we ask about 'true' platoons, 
defined as platoons of n > 2 vehicles. Using fl2.42p we immediately find that 
the density ptrpi(^^) of true platoons of velocity v is equal to 



Ptrpl(^^) = Plcad(^^)[l -Pl(^)] 

0-foll(^^) 



[p{v) - p{oll{v)] 



exp 



(2.43) 



Piv) - Pfoll(f ) 

We will not pursue these matters any further here. 
2.3 Circular traffic lane 

For the purpose of simulation or for theoretical considerations it may be more 
convenient to deal with periodic boundary conditions (PBC). We consider 
in this section the changes that this brings about with respect to the earlier 
case of open-ended boundary conditions (OBC). 

We consider one-lane traffic on a circular lane of circumference L. The 
spatial density p{v)dv = pR{v)dv of vehicles with natural velocity in [v,v + 
dv] is now the prescribed control function. We write again w = 0(f) for the 
effective velocity of a vehicle of natural velocity w. In this case the question 
is to find an equation that expresses 0(f) in terms of p{v). A shortcut way to 
do so is to assume that imposing boundary conditions along a different axis 
does not lead to a new type of stationary state and that the key equation 

pW = ^. (2.44) 

0(f) 

derived for OBC, remains valid for PBC with the same function 0(f). When 
we use fl2.44p to eliminate w(f ) in favor of p{v) from equation (12.181) we obtain 

(j){v) = v- I df'p(f')r(f')(f -f') [0(f) -0(^^')], (2.45) 

where p{v) is given and 0(f) is the function to be solved. In appendix [El we 
present a first-principle derivation of fl2.45p that does not make use of this 
shortcut. 



^This is most likely an approximation. 
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In order to solve equation fl2.45p we set 

= v+ dv' p{v')t{v'){v - v')(t){v'), 



vo 

V 



h{v) = 1+ / dv' p{v')t{v'){v-v'). 

J VQ 



(2.46) 
(2.47) 



in which <^iv) is unknown. Then 



^'{v) = l + j^dv' p{v')r{v')cP{v'), 
h'{v) = j:dv'p{v')T{v% 



<l>"(f) = p(f)r(f)0(f ), 
h"{v)=piv)Tiv), 



(2.48) 



which are such that /i(fo) = 1 and h'ivo) = 0. In terms of $ and h equation 
(I2.45P takes the form 

^ = m (249) 
^v) h{v) ■ ^ ' 

The desired solution must satisfy the boundary conditions $(fo) = vq and 
^'{vq) = 1. One obtains 

dv' 



$(f ) = h{v) 



Vo + 



(2.50) 



-.0 h'iV) 

as may be verified by direct substitution of fl2.50p in f l2.49p . Using f l2.50p and 
( ]2.48p we then arrive at the desired expression for 0(f) in terms of p{v), 



4>{v) 



£>) 
h"{v) 



Vo + 



dv 



h^{v') 



Vo + dv' 



vo 



1 + 



dv" p{v")T{v")iv' - V") 



. (2.51) 



All relations found before involving piead and pfon remain valid with uj{v) 
replaced by p{v)(f){v). In particular, one may rewrite fl2.32p as 



Plead (^^) = p(f ) 



dv' [(Piv) - (Piv')]p{v')T{v') 



VQ 



(2.52) 



By substituting fl2.5ip in f l2.52p one may eliminate 0(f) and obtains after 
simplifying 

p(f) 

This is the PBC analog of the OBC equation fl2.35p . 
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2.4 The moments X and Y of /9iead(^) 

The following notation will be useful later on. We will let X and Y denote 
the zeroth and first moment, respectively, of piead('^^), that is, 

X = dwpiead(^^) = Plead , (2.54a) 
Jvo 

POO 

Y= dvvp\cad{v) = P\ci,d{v)lcad, (2.54b) 

SO that (f)iead = Y / X IS the average platoon velocity. 

2.5 Upper density limit 

The model of this section obviously contains simplifications that may all be 
a matter of discussion. One of these, however, needs to be addressed in view 
of our ultimate conclusions. The fact that in this model the vehicles and 
platoons are all considered to have zero length, sets a limit on its validity 
when the total vehicle density increases. We may formulate the limit of 
validity as the condition that the typical platoon length ipi be much smaller 
than the typical distance (iinterpi between two successive platoons. Now we 
have ciinterpi = 1/piead, and may write £pi = ndintcrveh, where ciinterveh is the 
typical distance between two successive vehicles in a platoon (and includes 
the typical length of a vehicle). Hence the condition ipi <^ (iinterpi becomes 

^ ciinterveh ^ ~ • (2.55) 

Plead 

When with the aid of ( I2.39p we eliminate n from (I2.55p we obtain the limit 
of validity 

P ciinterveh < 1- (2.56) 

Here (imterveh is a parameter that is external to the present model (where it 
has in fact been set equal to zero), but may be estimated from observations 
of real-life traffic. We will check the status of this condition numerically in 
section 14. 4[ 

3 A model for two opposite lanes 

We now consider two lanes A and B with opposite traffic flows. In the 
general case of nonidentical traffic conditions on the two lanes there will be 
two control functions, UAiv) = UAPAiv) and u!b{v) = u}bPb{v), and two 
queuing time functions, TAiy) and tb{v). In the one-lane problem of the 
preceding section the queuing time function was considered given. In this 
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section we will aim at expressing the queuing time ta{v) experienced by the 
traffic in lane A in terms of the traffic conditions in the opposite lane B, 
and vice versa. First, in subsection I3.H we introduce the concept of the 
"encounter rate" between a vehicle in one lane and the platoons coming in 
the opposite direction. This builds straightforwardly on the work of section 
[2J Secondly, in subsection 13. 2[ we need a mean-field type hypothesis: we 
assume that when attempting an overtaking maneuver, each vehicle in one 
lane is subject to the average effect of all vehicles in the opposite lane. This 
will result in a set of equations in which the two lanes appear coupled. The 
structure of these equations is briefly discussed in section 13.31 From these 
coupled equations we will then in section H] be able to solve self-consistently 
- albeit only numerically - the functions ta{v) and tb{v). 



3.1 The encounter rates 7^ and 75 

As a preliminary let us consider two independent single-lane models A and 
B, one for each of the two directions. Let 'jAiv) [or ■Jb{v)] denote the rate 
at which a platoon of velocity v in lane A [or in lane B] encounters platoons 
in the opposite directioijfl; 'jAiv) and 7B(f) are necessarily an increasing 
function of v. The same methods as used in the single-lane model lead 
straightforwardly to expressions for these encounter rates, as we will see 
now. 

During a time interval [0, T] a marked platoon in lane A initially supposed 
in the origin will encounter all lane B platoons of velocity Vb that are initially 
in the space interval [0, {v+Vb)T]. Let NB,iead{v, Vh)dvb be the average number 
of platoons on lane B encountered by the marked one on lane A that have 
their velocity in [vb,Vb + dv^]. Then 

A'lead(f , Vb)dVb = iv + Vb)T X p B yi.d{vb)dvb . (3.1) 
The total number 7^(f )T of platoons encountered by the marked platoon is 

/■oo 

'yj^{v)T = dvbNB,ies.d{v,Vb) (3.2) 
and therefore from fl3.ip and f l3.2p 

POO 

1a{v) = dVb{v + Vb) pB,lei,d{Vb)- (3.3) 



Hence, in the notation of equation f l2.54p augmented with the appropriate 
lane indices, 

7a(i^) = PB,iead [v + {v)B,\ci,d] = XbV + Yg, (3.4a) 



■^We take all velocities positive, so that two vehicles of velocities Va and vt, in opposite 
directions have a relative velocity Ua + Wb- 
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and, by a permutation of the indices, 

1b{v) = PA.lead [v + (v) A,les.d\ = XaV + Ya- (3.4b) 

We conclude this subsection by a remark. 

Remark. - Let ^a,b be defined as the averages of the ^a,b{v) of equations 
(13.41) with respect to all platoons in lanes A and 5, respectively. Explicitly, 

lA = PB,lcad[(^^)yl,lead + (^)B,lead] = XbYa/ Xa + Yb, (3.5a) 
7B = PA,lead[(^^)A,lead + (f)i3,lcad] = XaYb/Xa + Ya- (3.5b) 

Then the expression 

PA,lead7A = P_B,lead7-B = PA,leadPB,lead [("y) A.lead + {v)B,lca.d] 

= XaYb + XbYa, (3.6) 

which is symmetric under exchange of the two lanes, represents the number 
of encounters per unit of time and unit of road length of platoons traveling 
in opposite directions. 

3.2 Coupling the two lanes 

In order to establish the equations that couple the two lanes, we reason as 
follows. Let To be the elementary time interval spent by an overtaking vehicle 
in the lane that is not its own. We consider tq as a fixed constant, in practice 
of the order of 15 seconds. We now consider that a vehicle queuing behind 
another one will overtake as soon as the platoons arriving in the opposite 
direction allow for a time interval at least equal to tq. The times ta{v) and 
tb{v) are the queuing times that result from this condition. Hence they are 
determined by tq and the encounter rates with the traffic in the opposite 
lane. We will now construct explicit expressions for ta,b{v). 

In order to relate 'yA,B{v) to ta,b{v) we introduce the additional hypothe- 
sis (briefiy commented upon in section [5]) that the platoons arrive at random 
times. The result, obtained in appendix |Fl then is that there is a function 
F{x), increasing monotonously from to oo on the positive real axis, such 
that 

^A,Biv) = roF{'yA,B{v)To), (3.7) 

where we recall that the 7^,5 are given by equation (13.41) . Appendix [F] yields 
the explicit expression F{x) = {e^ — l—x)/x, but in what follows no doubt any 
other F{x) with qualitatively the same behavior would lead to qualitatively 
the same results. 

Equation (13.71) is the desired relation that expresses the queuing times in 
one lane as a function of the rate of encounters with platoons in the opposite 
lane. It represents a coupling of the mean-field type: the traffic in one lane 
reacts to the average traffic conditions in the other lane. 
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3.3 The two- lane equations 



Let the lan e B a verages Xb = Pij,iead and Yb = PB,iead(w)s,iead be known. 
Because of f l3.4ap and (13. 7p we then know ta{v) in terms of Xb and Yb . We 
next use this expression for Tj^{v) in the single-lane theory of section [2J When 
substituted in equation fl2.35p [in the case of OBC] or in equation fl2.53p [in 
the case of PBC] it yields pA,iea.d{v), from which follow by f l2.54p the lane A 
averages Xa = p^.iead and Ya = PAAcad{v)AAead, again in terms of Xb and Yb. 
We may write concisely 

Xa = JaIXb, Yb), Ya = QaIXb, Yb). (3.8a) 

A permutation of indices yields 

Xb = fsiXA, Ya), Yb = QbIXa, Ya). (3.8b) 

This system of equations may be solved numerically, which will be done for 
a few examples in section m 

The indices on / and g are necessary because these functions depend on 
loa and cob, respectively, in the case of OBC, and on pa and pb, respectively, 
in the case of PBC. For equal control functions, UAiv) = ub^v) or Pa{v) = 
Pb{v) (which we will call the 'symmetric' case) we have fA = fB = f and 
9a = 93 = g and hence 

Xa = f{XB, Yb), Ya = g{XB, Yb), (3.9a) 

XB = f{XA,YA), YB = g{XA,YA). (3.9b) 

The symmetric case will be the one of main interest in what follows. If, in the 
symmetric case, we wish to just obtain a symmetric solution Xa = Xb = X 
and Ya = Yb = Y oi (13. 9p . we can simplify that equation further and get 
X = f{X,Y) and Y = g{X,Y). However, under certain conditions the 
symmetric case will be seen to allow for solutions with spontaneously broken 
symmetry; to find these the full equations (13. 9p are necessary. 



4 Numerical solution of the two-lane equa- 
tions 

The two-lane equations may be solved numerically. The procedure runs as fol- 
lows. We assume suitably chosen initial values for the pair {Xb, Yb) and then 
alternatingly update {Xa, Ya) and {Xb, Yb) by applying equations (I3.8ap and 
f l3.8bp . In practice equations (I2.35P and (I2.53p . which are at the heart of the 



iteration procedure, have to be discretized. We use the natural discretiza- 
tion, that is, the form that these equations take when the control functions 
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<^a,b{v) or pa,b{v) are sums of Dirac delta's. This discretized version is de- 
scribed in appendix [Cl where equation f l2.35p has as its equivalent in (1C.7P 
together with (]C.5|) . and f l2.53p has as its equivalent in (1C.13P together with 

In all cases that we have considered there is convergence to a solution, 
which appears either as a fixed point with (X^, Y^) = {Xb, Yb) or as a fixed 
two-cycle alternating between {Xa, Ya) and {Xb, Yb). The convergence slows 
down near the critical point (see below), but we have not found it worthwhile 
to improve the algorithm. 



4.1 Symmetric lanes with OBC 

We consider first the symmetric case with identical a priori traffic conditions 
for the two lanes, that is, equal distributions uja{v) = u)b{v), which we may 
therefore denote by (xj{v) = ujP{v). We recall that uj is the total intensity of 
the traffic fiow, that is, the rate of entrance of vehicles, irrespective of their 
velocity f , at the point x = of the lane section under consideration. For the 
distribution uj{y) we take a truncated Gaussian discretized on N + 1 points 
Vi = Vq + iAv with i = 0,1, . . . , N. Explicitly, 



AT 



^UiS{v-Vi) (4.1a) 



with oji = upi and 



Pi 



ATexp (^-^^^i^-^) , ^ = 0,1,...,N, (4.1b) 



where Vc is the peak velocity and Af the normalization constant required to 



satisfy J^iLoPi = 1- We have studied numerically the example with 
fo = 60km/h, cr=10km/h, Af = 1 km/h, 

(4.1c) 

t;c = 90km/h, N = QO, 

so that the largest occurring velocity is fm = 120km/h. When averaged over 
all vehicles in a road section, and in the absence of interaction between the 
lanes, equation (14. Id) yields a mean vehicle velocity of (v)^'^^ = 88.9 km/h (it 
is less than the average velocity of 90 km/h since slower vehicles stay longer 
in the section than faster ones). 

For this example we determined the solutions {Xa,Ya) and {Xb,Yb) 
with the control parameter u varying from w = to about u = 1500 ve- 
hicles/hour. Figure [7] shows the results represented in the plane of coor- 
dinates {X,YX~^) = (piead, ('J^)iead)- In the branch marked 'A,B' the pa- 
rameter Cj runs from zero in the upper lefthand corner to a critical value 
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uj = ujc = 800.5 h^^ at a point where the curve bifurcates. The interpretation 
is that for u < Uc the two lanes A and B have identical behavior; but when 
UJ reaches its critical value, the symmetric solution becomes unstable and the 
symmetry between the lanes A and B is spontaneously broken. From the bi- 
furcation point two branches marked A and B emanate; in the figure Cj runs 
from Uc to a value of the order of 1500 vehicles/hour along these branches. 
Lane B has a high total density of vehicles, many of which are included in 
platoons, and the other lane A has a lower density of vehicles, the majority of 
which are leading, and thus many of them must constitute 1-platoons. The 
continuation of the 'A,B' branch, where the symmetric solution is unstable, 
has been indicated dashed line. 

Figure [8] shows, as a function of the control parameter u, the densities 
Plead and pfou of the leading and following vehicles, as well as their sum p. 
Figures [10] and [11] show the corresponding graphs of the average platoon 
length n and the average platoon velocity (v)iead- For small u we are in the 
weakly interacting regime. The three densities shown in figure [S] are linear 
in cu, with deviations from linearity becoming appreciable when u attains 
values of a few hundred vehicles/hour. We recall now that piead is also the 
density of platoons. Figure [8] shows that for small u the platoon density piead 
is practically equal to the total vehicles density p, which means that most 
vehicles constitute 1-platoons by themselves. This is confirmed by the initial 
behavior of n in figure [10] Consistently, (f )iead in figure [11] stays initially 
equal to its noninteracting value {v)^'^\ 

Beyond the small u regime the average platoon velocity begins to decrease 
due to faster vehicles getting trapped in slow platoons. The instability occurs, 
in this example, when the average platoon length has increased to about 
two cars per platoon and the average platoon speed has gone down to about 
85 km/h. The interpretation of this symmetry breaking is most easily derived 
from the large Co behavior of the two branches. As u increases to values well 
above Uc, there occurs formation of a few very long platoons behind vehicles 
having velocities close to the minimum velocity. This explains why on branch 
B asymptotically piead and (f )iead Vq. The small number of platoons 
in this lane have large spaces between them which allow the vehicles in lane A 
to overtake almost without having to queue. Hence asymptotically on the A 
branch (f )icad ~^ {v)^'^\ This interpretation is again confirmed by the figures 
[5|[TUl and[TTJ Particularly striking is the rapid increase of the average platoon 
length in lane B when u increases in the supercritical regime. It returns to 
the average vehicle velocity in the absence of interaction; in lane B it tends 
asymptotically to the velocity vq = 60 kmh~^ of the slowest vehicles. 

We note that the bifurcation points of piead and pfou seem to coincide in 
figure [HI but figure [3 shows that they are actually distinct. 

Finally in figure[T2]we show the behavior of the effective velocity w = (j){v) 
for two values of u. They all have the expected behavior. For the subcritical 
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Figure 7: (OBC) For the velocity distribution u}{v) = CjP{v) given by equations 
(14. ip . as the control parameter a) is varied, the two-lane equations (I3.9P yield a 
locus of solutions of solutions has been plotted here in the plane of coordinates 
{X,YX~^) = (piead, (^^)iead)- The control parameter increases along the branch 
marked 'A,B' from a) = in the upper left corner to around a) = a)c ^ 800 h-i at 
the bifurcation point, and, going outward from this point, from a) = Wc to around 
u) = 1400 h~^ along each of the two branches marked A and B. The dashed line 
represents schematically a branch of unstable solutions. 



flux Oj = 500 vehicles/hour the functions (pAi^) and (psiv) coincide. For the 
supercritical value Cj = 1000 vehicles/hour they are distinct; the lane A curve 
is not far from the straight line 0(f) = v valid for the noninteracting case, 
and in the lane B curve the effective velocities are reduced to values not far 
from the minimum velocity Vq = 60km/h. 

We have varied the distribution P{v) and found that this does not affect 
the qualitative features of the results that we discussed. 

4.2 Symmetric lanes with PBC 

It is worthwhile also to see what happens when the open boundary conditions 
(OBC) are replaced by periodic boundary conditions (PBC). For one thing, 
PBC may be easier in simulations of this and other models. In the case of 
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Figure 8: (OBC) The total vehicle density (marked T), which is the sum of the 
leading vehicle density (L) and the following vehicle density (F), plotted as a 
function of the flux oj for the same distribution uj{v) (common to both lanes) as 
in figure [71 At the critical flux a)c = 800.5 vehicles/hour each of these three curves 
bifurcates into an A and a B branch, signaling spontaneous symmetry breaking 
between the two lanes. 

PBC the density p{y) = pR{y) is given and the total vehicle density p is the 
control parameter. We have studied the example where R{v) is the same 
truncated Gaussian, discretized on + 1 points f j, that was used in section 

sn 

The flux u is now a derived quantity. We may decompose it as the sum 
of the fluxes cjiead and cjfou of leading and following vehicles, respectively. In 
figure US] we have plotted u and cjiead as functions of p. Both bifurcate at a 
critical point p^ = 9.02 km~^. 

The behavior of the densities piead and pfou may also be found as a function 
of p; it is very similar to the case of OBC, but with the difference that their 
sum p is now a nonbifurcating straight line of unit slope. It thus appears 
that OBC and PBC lead to the same qualitative pictures. The qualitative 
features again remain unchanged as one varies the distribution R{v). 
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Figure 9: Detail of figure [HI Closed squares: /Oiead- Open circles: /Ofou. 
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Figure 10: (OBC) Average platoon length as a function of the flux a) for the same 
distribution u}{v) (common to both lanes) as in figure [71 For O > uic the platoon 
length in lane B is seen to increase very rapidly with the flux. 
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Figure 11: (OBC) The average platoon velocity as a function of the flux u) for 
the same distribution oj{v) (common to both lanes) as in figure [71 For u) > lo^ the 
platoon velocity in lane B is seen to decrease rapidly with the flux. 
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Figure 12: (OBC) The effective velocity (f){v) of a vehicle as a function of its 
natural velocity v. The distribution uj{v) is the same as in figure [7] and common 
to both lanes. At the subcritical flux of a) = 500 vehicles/hour lanes A and B 
are described by a common curve, whereas at the supercritical flux Co = 1000 
vehicles/hour they are described by two separate curves. 
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Figure 13: Case of periodic boundary conditions (PBC) with vehicle density 
p{v) = pR{v) and R{v) a truncated Gaussian described in the text. Shown are the 
total vehicle flux u) (marked T) and the leading vehicle flux wiead (marked L) as a 
function of the vehicle density p, which is here the control parameter. 
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Figure 14: (OBC) Locus of solutions of equations (j3.9p in the asymmetric case 
with ujb = 0.95ti)A- Branches A and B go over continuously into those marked Al 
and Bl, respectively. At a critical value Qa = ^A,c ~ 901 vehicles/hour a secondary 
solution appears with two stable branches marked A2 and B2. This solution has 
B as the fast lane and A as the slow one. The dashed lines schematically represent 
unstable parts of the secondary branches. 
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Figure 15: The leading veliicle density piead corresponding to figure [TH At ua = 
WA,c = 901 vehicles/hour two new branches appear, marked here LA2 and LB2. 
The unstable branches are no longer shown. 
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Figure 16: Average platoon length n corresponding to figure [TH The unstable 
branches are not shown. 
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4.3 Non-symmetric lanes with OBC 

There are many ways to break the symmetry between the two lanes. We have 
chosen to investigate the open-ended problem with ua 7^ ojb, Sill other things 
in the two lanes being identical. In particular, Pa{v) = -Pb(v) = Piv) is again 
the Gaussian distribution of equations (]4.1b|) - fl4.1cp . We have taken a small 
symmetry breaking represented by ub = 0.95 cua and solved again equations 
f l2.43p for Xa,Ya, Xb, and Yb- Instead of figure [7] we now obtain figure [HI 
in which again the dashed lines represent (only schematically) branches of 
unstable solutions. Correspondingly, instead of the curves marked L, LA, 
and LB of figure E] we now obtain the curves in figure [T51 For low values 
of ua (and concomitantly ub) the platoon densities in the two lanes (the 
curves marked LA and LB) increase linearly with ua, their slopes having a 
ratio 1 : 0.95. When Cja increases, the critical point is avoided and the platoon 
density of lane A (lane B) evolves continuously to the values represented by 
the branch marked LAI (LBl). However, there is now a new critical value of 
uja = ^A,c =901 h~^ at which an inversion of the densities becomes possible. 
The two lanes are then on the 'secondary' branches LA2 and LB2. 

We are in the presence of a stationary state problem for which there is 
no obvious way to define a free energy functional, although this may not be 
impossible. However, it seems reasonable to consider this set of secondary 
branches as 'metastable' (in the sense of having a higher free energy) with 
respect to the main branches. 

The curves for the total and the following vehicle densities (not plotted 
in figure [T^ have analogous behavior. In figure we have plotted the 
corresponding average platoon lengths. 

4.4 Upper density limit 

As we have seen in section 12. 5[ there is a high density regime where the 
present model ceases to be valid. The condition for validity depends on the 
external parameter ciinterveh, which represents the typical distance between 
two consecutive vehicles in a platoon. For the example of equations (14. ip we 
may estimate rfinterveh = 25 m, which when used in equation fl2.56p gives 

p < 40 vehicles/km. (4.2) 

As we have seen, the principal phenomenon investigated in this work, namely 
the spontaneous symmetry breaking between two traffic lanes, occurs at a 
density of p ~ 5 vehicles/km. Our conclusions concerning this phenomenon 
are therefore compatible with the simplifications of zero vehicle length and 
zero platoon length that we used. 
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5 Mathematical remarks 



In this penultimate section we collect some comments on the calculations 
presented at various stages of this work with the purpose of clarifying their 
precise mathematical status. 

The one-lane model is a statistical ensemble of geometrical constructs 
of the type represented in figure [61 This ensemble is defined by (i) the 
construction rules of sections [1] and [21 and (ii) as an input, the statistics of 
the intersections of the vehicle trajectories with one of the boundaries of the 
space-time domain, either the t or the x axis. One of the simplest ways of 
realizing (ii) would be to assume Poisson statistics for the intersection points 
and to assign to the vehicles velocities drawn independently from a given 
distribution, whether P{v) or R{v). More generally, however, one could on 
the boundary specify arbitrary correlations between the positions and the 
velocities of the vehicles. 

To develop the theory we have imposed, in addition to P{v) or R{v), only 
the average flux Cu or average density p, but have not specified any boundary 
correlations. Our theory is therefore essentially a set of relations between 
average densities, fluxes, velocities, etc. 

In section 12.1.31 the existence of a stationary state and a well-defined 
function 0(f) is strictly speaking an assumption, however natural it may be. 
We leave open the possibility that the stationary state correlations (but not 
the averages) depend on the unspecified boundary correlations. 

On only two occasions did the stationary state correlations between the 
vehicles play a role and was an extra assumption about them needed. The 
first occasion was when we discussed the statistics of platoon lengths in 
section I2.2.4( this calculation was accessory and had no bearing on the rest 
of the paper. The second occasion was in section 13.21 and appendix [Fl when 
we derived the function F that connects the queuing times in one lane to 
the traffic density in the other one; although this connection is the very 
essence of the two-lane model, we expect that different assumptions about the 
correlations would lead to a qualitatively similar function F. Nevertheless, 
in the single-lane model and for fully specified boundary statistics, it is a 
well-defined and still open question to determine the statistics of the platoon 
arrival times at a fixed observation point. We have not attempted in this 
work to answer that question and have simply presented the Poisson process 
nature of the arrival times as an additional hypothesis. 

6 Summary and conclusion 

In this paper we have proceeded in two steps. 

First we defined a geometric traffic model of fiux ua on a single lane 
('lane A'). The model includes the possibility for fast vehicles to overtake 



38 



slower ones and is attractive because of its simplicity. Vehicle trajectories are 
straight lines except that overtaking vehicles incur time delays that follow 
from a given 'queuing function' ta{v')^ where v' is the velocity of the vehicle 
being overtaken. The single-lane model exhibits formation of platoons of 
vehicles. The platoon density piead as compared to the total vehicle density 
p has been one of our main quantities of interest. We were able to relate 
analytically the platoon's statistics to the flux intensity. 

Secondly, and on the basis of additional hypotheses, we expressed the 
queuing function ta governing the lane A vehicles in terms of the traffic con- 
ditions in the opposite lane, 5, and vice versa. This resulted in a set of cou- 
pled two-lane equations. Although the two lanes A and B have individually 
been modeled microscopically - that is, in terms of the trajectories of individ- 
ual vehicles -, the two-lane equations connect densities and average velocities 
in the two lanes but no longer have a microscopic interpretation. We solved 
the two-lane equations numerically for a characteristic example of traffic flow 
with a priori identical traffic conditions on the two lanes {uja = ojb = It 
was found that above a critical traffic flux, uj > COc, (or above a critical traffic 
density) the symmetry between the two lanes is spontaneously broken. One 
lane has fast traffic with mostly 1-platoons (= single vehicles), whereas the 
other lane has slow traffic with almost all vehicles queuing in long platoons. 
Asymmetric conditions {ua 7^ ^b) have been briefly investigated, as well as 
the case where not the fluxes oja,b, but the spatial vehicle densities pa,b are 
the control parameters. 

We argued that the model considered in this work is valid for low and 
moderate traffic densities. A more elaborate version, appropriately defined 
and including the parameter d^^terveh (= the typical intervehicle distance in 
a platoon), would be able to describe also the very high density regime and 
the phenomenon of jamming of one or both lanes. However, the spontaneous 
symmetry breaking reported in this paper occurs in the range of densities for 
which the present model is valid. 

It has not been our purpose to be exhaustive. Several questions of interest 
have been left aside here but can be fairly straightforwardly treated in the 
same model. One of these concerns cases in which the velocity distribution 
has a tail that goes down to zero. Another one is the incorporation, along the 
road section under investigation, of entrance or exit ramps with prescribed 
vehicle fluxes. A less simple extension would be to include time-varying 
boundary conditions as may be due, for example, to traffic lights at the 
entrance of the road section. We have considered only the stationary state; 
studying the relaxation of an initial state towards stationarity would allow, 
among other things, a comparison with the scaling laws of reference |32] . 

Analytically various limit cases could be pursued further. In the weakly 
interacting limit an expansion of all relationships of this work in powers of 
uta and utb is certainly possible. In the limit where the variance cr^ of 
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the velocity distribution becomes small, the phase transition point moves to 
higher densities; an expansion for small variance seems more complicated 
but nevertheless also possible. We have not exhibited any exact solutions to 
our equations; no doubt a few exist and we think some of these may be of 
interest. 

The lane-lane coupling that we have introduced in the present study is of 
mean-field type, since each lane A vehicle feels only the average effect of all 
lane B vehicles, and vice versa. One should wonder if this approximation is 
justified and if in its absence the symmetry breaking still exists. We believe 
it does, the main argument being that each vehicle in one lane encounters in 
the course of time all vehicles in the other lane. This argument is admittedly 
heuristic. Some light is shed on its validity in two ways. First, in sepa- 
rate work [35] we found that the simulation of a closely related microscopic 
two-lane model also presents this symmetry breaking, at least on practically 
relevant scales of time and space. Secondly, in the quite different context of 
magnetic friction Hucht and coworkers [311 EZ] were interested in an Ising 
system on a ladder lattice whose two legs are in relative motion. This model 
has essentially the same features as the present road traffic model; an exact 
solution [37] shows that there is a phase transition with mean-field charac- 
teristics for infinite relative velocity, which however gets rounded when the 
velocity becomes finite. We will leave further discussion of these issues to 
the future. 



A Solution of Eq.[239] for 

In this appendix we solve equation (12.191) for 0(f). We denote the numerator 
of its RHS by 



Upon setting 0(t>) = l/x{'^^) we find that x(^) satisfies the linear equation 




(A.l) 



This quantity satisfies 





(A.2) 




(A.3) 



whence by differentiating twice we obtain 




(A.4) 
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and 

[xg]" = 9"X. (A.5) 
Setting now ip{v) = ) we can rewrite f lA.5|) as 



[log^4']' = -2[\ogg]'. (A.6) 
Two successive integrations applied to (IA.6P yield 

X{v) = ^{v) = (A.7) 

and 

-i- = xiv) = C, + C, f\v' , (A.8) 

where Ci and C2 are constants of integration. In order to determine these 
constants we first observe that (lA.Sp leads to l/0(fo) = x('^o) = C2. Since 
for the reasons exposed in section [T] we have 0(fo) = Vq, it follows that 
C2 = l/fQ. Next we evaluate (IA.4P for v = vq. Using that g{vo) = Vq and 
9'i'Vo) = 1 this leads us to foX'(^o) = ^x(^o) whence x'(^o) = ~l/'^o- Upon 
comparing this last result to (IA.7P evaluated at f = f we see that Ci = — 1. 
Substitution of the values of Ci and C2 in flA.Sp yields the solution f l2.20p 
given in the main text. It is easily verified that for u}{v) = equation f l2.2Up 
reduces to 0(f) = v, as had to be the case. 



B An analytic example of functions (j){v) and 

Aead(^; 



B.l Analytic example of (l){v) 

Let uj{v) = u}P{v) be a block distribution of the natural velocities, 

{OJ 
Vo<V < v^, 
else, 

let t(v) = f, and define the dimensionless parameter A = uf. Integrating 
according to equation (lA.ip we first find 

,,(„) = „+Ai^i^. (B.l) 

To abbreviate the notation we set 

k{u) = Vq + U + CV?, U = V — Vq , 

k'{u) = l + 2cu, c= X/{vm-Vo) , (B.2) 

A = -1 + 2foA/ (fm - Vq) . 
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Next we find from equation (]A.8|) tliat 

1 1 



(f){v) Vo 

wfiere I{u) is tlie indefinite integral 

du 



+ I{v-Vo)-I{0), 



I u) 



k\u) 



k'iu) 4c 
1- < 



1 k'iu) 
arccot- 



1 



Ak{u) A 



1 



artanh 



arcoth 



'A 

k'ju) 

k'ju) 



(B.3) 



A > 0, 



A<0, ^m>l. 



(B.4) 



Depending on the values of the parameters and the variables of the problem, 
all three cases of this equation may occur. The desired function (piv) now 
follows by substitution of (IB.4P in fIB.Sp . 



B.2 Limit of small A 

The above result becomes easier to visualize in the limit of small A. In that 
limit the 'arcoth' applies in equation flB.4p . After a certain amount of algebra 
one obtains 



[v) = V + 



A 



(3w — vq){v — vq) — 2f ^ lo^ 



^^0 



+ 0(A2), (B.5) 



in which the expression in brackets may be shown to be negative for all 
V > Vq, as had to be the case because of 0(f) < v. Expanding the coefficient 
of the term linear in A near the lower limit Vq of the velocity distribution we 
get 



(f){vo + u) = Vq + U- X 



u" 



+ 0{u' 



+ 0{X^). 



(B.6) 



Let {w)\ = (0(f))A denote the effective velocity w averaged with repect 
to P{v) at a given A = uf. Then we find from (IB.ip together with (]B.6P 



{w)x = {v)o + 



A 



2(fm - ^^o) 



if^m - 2^>0 + VmV^Q - - log 



+OiX^) 
(B.7) 
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in which {w)o = {v)o = |(t'm — vq) is just the average natural velocity. In the 
limit of a narrow distribution, f m — — ^ 0, this becomes 

{w)x={v)o - [iv,^-voY + o(^iv^-vof)]+0{\'). (B.8) 

The factor {v^—vqY shows that the reduction of the average effective velocity 
is proportional to the variance of the distribution. 

B.3 Analytic example of piea.d{v) 

The platoon density piead{v) corresponding to the present example is easi- 
est to obtain from equation f l2.35p [which is also (ID.6P -( lDl2|) ]. from which 
(f){v) has been eliminated. Upon observing that K{v) = k{v — vq) we get 
straightforwardly from (1D.6P the simple expression 

from which many further analytic results may be derived by means of the 
various relations of this work. 



C Discretized equations 

One may discretize the velocity axis by setting Vi = Vq + iAv, where Av is a 
fixed increment. We let i run from to some large value that may later 
on be set equal to infinity. The discretized version of the equations is used 
for the numerical solutions discussed in this work. It is, however, more than 
just a way of approximating the continuum limit, since it describes a traffic 
model having only a discrete set of velocities. 

Open boundary conditions. Let us now take for uj{v) = ujP{y) the discrete 
distribution 

N N 

^{v) = ^Ui5{v -Vi), ^Ui = u, (C.l) 

j=0 j=0 

and set 

r. = r(v,). (C.2) 
Upon substituting ( IC.ll) and ( 1C.2I) in ( 12.19^ and abbreviating 0j = (pivi) we 



get for 0j the equation 



' l + ^YU{j-i)ujfr,ct>r^ 



j = l,...,iV, (C.3) 



5 
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in which Tj = T{vi) and where the right hand member depends only on 
00, • • • , 4>j-i- Given that 0o = "^o we may obtain from this equation succes- 
sively the numerical values of 0i, 02, • • • , 0Af- One can easily see that they are 
increasing with j. Equation (IC.3P is the discrete counterpart of fl2.20p . Its 
explicit solution 0^ = 0(fj) is 

0,- = At; V (C.4) 

with 

i-i 

Qj = Vj + Av ^(j - i)uJiTi . (C.5) 

1=0 

When the velocities have the discrete distribution fIC.ip we may set 

N 
j=0 

and find that the discrete equivalent of (12.351) is 

Plead J = (^j9j^, j = 0, 1, . . . , A^, (C.7) 

where we use the notation of appendix O 

Periodic boundary conditions. In the case of cyclic boundary conditions 
the density p{v) is the control function. Let us suppose it is of the discrete 
type 

N 

p{v) = Y,Pr^{w,). (C.8) 

i=0 

Substitution of this expression in (I2.52p yields again expression fIC.Sp but with 
ujpi replaced with pj0j. Given that 0o = vq we may determine successively 
01, . . . , 07V in the same way as for open boundary conditions. 
Solving the equation explicitly yields 

^ 1 

0, = t;o + At;V— — , j = 0,l,...,iV (C.9) 

1=1 

with 

h, = l + Av - ^)P^r^ , J = 0, 1, . . . , iV. (C.IO) 

i=0 

These two equations are the discrete versions of (I2.5ip and (12.470 . respec- 
tively. 
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In the case of cyclic boundary conditions the expression for the density 
p(f , v') is again given by (I2.27P but with uj replaced with p. Its discrete 
equivalent, with an obvious definition of p^, is 



Pij = PiPjTj[(pi 



(C.ll) 



Similarly, equation fl2.32p has, in the case of cyclic boundary conditions, the 
discrete equivalent 



Plcadj' Pj 



i-1 



1 - ^pinK^j - 



i=0 



(C.12) 



When fIC.QP is used in flC.121) we obtain after substantial rewriting the sim- 
plified result 

Pleadj = Pjh'^^ (C.13) 

with hj given by fIC.lOl) . Equations f lC.lSP and (IC.lOP constitute the discrete 
equivalent of f l2.53p and f l2.47p combined. 



D Rewriting Eq. (ESI) 

We consider here expression fl2.32p for the platoon density piead. 



Plead (f) = (^{v) 



dv' 



uj{v')r{v')\ , (D.l) 



and show how it may be cast in a much simpler form. It is useful to abbreviate 

(D.2) 



Hence 



and, for v' < f , 



1 



/ dv'{ 


V — v')u{v')t 


J VQ 




1 


r dv. 






1 


r dv. 







Using ( 1D.4P in (ID.ip we obtain 

1 



Plead (t^) = (^(v) 



co{v')t{v') 



1 r 1 

~ dv' 



(j){v 



K^{v') 



df2 w(t;2)r(f2) 



(D.3) 



(D.4) 



(D.5) 



Vo 
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where to arrive at the second hne we performed an integration by parts. 
When with the aid of (]D.3P we ehminate the on the right hand side 

of (ID.SP we obtain 




1 dK{v'] 



dv' 



(D.6) 



where we used that K{vo) = Vq. Upon substituting for K{v) in the last 
hne of (ID. 61) its definition (ID. 21) one obtains the final result of this rewriting, 
equation ( I2.35P of the main text. 



E An equation for (f){v) in the case of PBC 

For the case of a single traffic lane with periodic boundary conditions, studied 
in section 12.31 we wish to find an expression for ) in terms of the control 
function p{v). The reasoning resembles the one of subsection 12.2.21 but is 
still sufficiently different that we repeat the main steps here. Let T{v) be the 
time needed for a vehicle of natural velocity v to cover a distance L, that is, 
to complete one full turn of the ring. Then 

n^) = -^ = -- (E.l) 

(p[v) w 

We consider a marked vehicle of natural velocity v and hence of effective 
velocity w = (p{v). Let v' < v and w' = (f>{v'), so that w' < w. Let 
N{w,w')dw' denote the average number of vehicles with velocities in the 
interval [w', w' + dw'] overtaken by the marked vehicle. To find an explicit 
expression for this number, we notice that the total number of vehicles on 
the ring having a velocity in the given w' interval is Lp{w')dw' . Of these, the 
fraction overtaken by the marked vehicle is equal to 1/L times the relative 
distance {w — w')T{v) covered, that is, 

N{w,w')dw' = - w')T{v) ^ ^p^^,^^^, 

Lj 



w — w' 



= L p(w')dw'. (E.2) 

w 

As before, the time lost in overtaking a vehicle of velocity v' will be denoted by 
r(f'). Then, due to its overtaking of vehicles with velocities in [w' ,w' + dw'], 
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the marked vehicle, while it cycles the ring once, will be queuing during a 
total time Tfo\\{w,w') that is given by 



Tfo\\{w,w')dw' = t{v')N {w ,w')(kv' 



(E.3) 



(in hybrid notation where we use both v and w variables). The total distance 
dio\\{Wi w')dw' that it covers during this lost time is 



d{oii{w,w')dw' = T{v')N{w,w')v'(kv' . 



(E.4) 



Integrating now over all w' we find that the total time loss Tfoii(w) suffered 
by the marked vehicle during one cycle is 



Tfou (w) 



dw' t{v')N{w, w'). 



(E.5) 



Similarly the total distance dio\\{w) covered during this lost time is 



dio\\{w) 



dw' t{v')N{w,w')v' . 



(E.6) 



In the remaining time interval T — Tfoii(^) the remaining distance L — dio\\{w) 
is covered at a velocity v. Hence we must have the key equation 

v[T- ffou(«;)] =L- d^iiw). (E.7) 

We substitute flE.5|l and flE.6|) in this equation, using for N{w, w') the explicit 
expression flE.2p . The result is the relation 

: 1 - 



1 

w 



dw t{v ) p[w ) 



VQ 



W 



dw'T{v')^^^^v'piw'). (E.8) 



Do 



W 



When we divide by v and transform from w and w' to v and v' using (12. 9p 
we get 



1 



(v) V 



dv'T{v') 1 



l--]p{v'). 



(E.9) 



We can rewrite this as (12.451) . which completes the proof and indirectly con- 
firms the validity of (12.441) for the case of PBC 



F Relation between TA,Biv) and 7a,s(^) 
F.l A problem in statistics 

The question of finding a relation between ■yA,B{v) and ta,b{v) leads us to 
the following problem in statistics. Point events take place at a rate 7. We 



47 



start observing them at an initial time t = 0. Let the stochastic variable f 
denote the instant of time at which the first time interval starts that has a 
length greater than tq and is free of events. The probability distribution p{f) 
of f will be of the form 

p{r) = ASir) + B{r), (F.l) 

where A = e~''"^° is the probability that there is no event in [0, Tq] and where 
B{f )df is the probability that the first eventless interval greater that Tq 
starts in [f,f + df]. The calculation of B{f) is possible but the expressions 
become more complicated than is needed for our purpose. We will therefore 
do something simpler, but good enough for our purpose, namely replace f 
by its average 

fOO 

fj9(f)df. (F.2) 





For dimensional reasons we can write 

r = roF(7ro). (F.3) 
Furthermore we should necessarily have 

F(0) = 0, F(oo) = oo. (F.4) 
It is not difficult to find an explicit expression for r as defined by flF.2l) . 



F.2 The expression for rA,B{v) 

Let pq be the probability that there is no event in [0,ro], and let p„ be the 
probability that the first time interval of length greater than tq comes along 
after the nth event, for n = 1, 2, . . .. Then 

Pn = {l- e-'^^r e-^"", n = 0, 1, 2, . . . (F.5) 
Using ( ]F.2|) we get 

oo 

E^^avPn (F.6) 



oo 

r 

n=0 



where is the average length of a time interval given that it is shorter than 
Tq. An elementary calculation gives 

r< = / dt^ (F.7) 

1 — e"''^'' — 'yToe'^^" 
= 7(l-e-7^o) 



48 



When (IF.SP is substituted in (\FM we get 



E 

n=0 



:,7f"0 



1) 

1 -7^o]/7, 



(F.9) 



which is (13. 7p . It is indeed of the form (lF.3P - (]F.4p with F{x) = {e^ — l~x)/x. 
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